AREA LAW FOR RANDOM GRAPH STATES 

BENOiT COLLINS, ION NECHITA, AND KAROL ZYCZKOWSKI 

Abstract. Random pure states of multi-partite quantum systems, associated with 
arbitrary graphs, are investigated. Each vertex of the graph represents a generic 
interaction between subsystems, described by a random unitary matrix distributed 
according to the Haar measure, while each edge of the graph represents a bi-partite, 
maximally entangled state. For any splitting of the graph into two parts we consider 
the corresponding partition of the quantum system and compute the average entropy 
of entanglement. First, in the special case where the partition does not "cross" any 
vertex of the graph, we show that the area law is satisfied exactly. In the general 
case, we show that the entropy of entanglement obeys an area law on average, this 
time with a correction term that depends on the topologies of the graph and of the 
partition. The results obtained are applied to the problem of distribution of quantum 
entanglement in a quantum network with prescribed topology. 
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1. Introduction 

Entanglement in many-body quantum systems is a subject of a considerable recent 
interest [THS]. In several physical problems it is important to describe correlations 
between two selected parts of a composed quantum system. For any pure state de- 
scribing the entire system, such correlations can be characterized quantitatively by the 
entanglement entropy H, equal to the von Neumann entropy of the reduced mixed 
state obtained by the partial trace over a selected subsystem. Note that this quantity 
depends explicitly on the partition of the entire system into subsystems. 

l 
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Entanglement entropy H does not usually scale proportionally to the size of the 
selected region S. As described in [3111], for different systems the entanglement entropy 
is approximately proportional to the boundary of this set, denoted by dS. Considering 
a three dimensional body S, the size of its boundary is proportional to the area of S 
and not to its volume. Thus, for any set S of subsystems, the size of its boundary \dS\ 
will be called the "area" of S. If for a sufficiently large system the leading contribution 
to the entanglement entropy of a given state |^) with respect to the partition {S, T} is 
proportional to the area separating both subsystems, we will say that the area law is 
satisfied. In the case of one dimensional systems, the area law implies that the entropy 
saturates asymptotically to a constant, as the area of the boundary consists of two 
isolated points. 

Area laws are studied in context of black hole physics and holographic principle J5JE] , 
and for ground states of quantum lattice systems with local interactions [7HTU] . Under 
some technical assumptions is it possible not only to show that the area law holds 
for the ground state of a given two-dimensional model system, but also to derive the 
negative correction term called topological entanglement entropy [TTUT2] , and show that 
it is universal and depends only on the topology of the interaction. Investigations of 
the area law were performed also for subsets with fractal boundaries [13] . 

For any quantum state of a system of interacting particles, described by a lattice or a 
graph, it is interesting to analyze the degree of entanglement between two given nodes 
of the graph [2] . This issue is relevant in studies on generating entanglement between 
given nodes of a graph [15], entanglement swapping and quantum repeater systems [T] 
and quantum communications in noisy networks |16j . 

The aim of this work is to study the area law for the ensembles of random states 
associated to a graph. The authors have previously introduced [17] an ensemble of pure 
quantum states for which the structure of interactions and entanglement are encoded 
in the graph. In the present work, we study the average entropy of entanglement for 
elements of this ensemble, and we show that is obeys asymptotically and on average 
area laws of the form 

(1) E{V) = H( Ps ) = \dS\ log N - h r>s + o(l), 

where \I> is a graph state, ps is the reduced density operator, \dS\ is the "area" of 
the boundary of the partition {S,T}, N is the dimension of the subsystems and hr,s 
is a numerical constant depending on the graph T and on the partition {S,T}. The 
functionals E and H denote, respectively, the entropy of entanglement for pure states 
and the von Neumann entropy of mixed density matrices. 

In previous papers one was usually studied entanglement entropy for a concrete quan- 
tum state of a given composite quantum system - see e.g. [3J. Here however, we discuss 
the statistical properties of ensembles of quantum states. We will show under which 
assumptions the area law holds exactly in our model, while in the opposite case the 
corrections to the area law are derived. 

The paper is organized as follows. The model of graph random states from [TTJ is 
recalled in Section [21 and a detailed example is worked out. In Section [3j the notion of 
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adapted marginals of graph states is introduced and it is shown that for such partitions 
of the graph the area law holds exactly. A more general case of the problem, for which 
the area law for random graph states holds only asymptotically is treated in Section HI 
In Section El a definition of the boundary is proposed in the most general setting. The 
general area law is stated in Section |6l while an application to a transport problem is 
presented in Section [71 Section [8] provides an outline of the adjustments to be made 
to the main results in the case where the dimensions of the relevant Hilbert spaces are 
different. 



2. Random graph states 

We recall here the model of random graph states introduced in [17] . The most general 
definition of this model will not be recalled here, but for completeness we sketch below 
the main ingredients of the construction. For any undirected graph T consisting of m 
edges Bi, . . . , B m and k vertices v\, . . . Vk, we associate an ensemble of random pure 
states They describe a quantum system consisting of n = 2m particles described 
in the composed Hilbert space H% 2m . The dimension N of any subspace is arbitrary 
and in particular we will analyze the asymptotic limit iV — > oo. 

Each edge of the graph, which connects subsystems labeled by i and j, represents 
the maximally entangled state, \$)fj = ^= XL=i between two Hilbert spaces 

Hi and %j. Each vertex V of the graph, of degree b, represents a generic interaction 
between b subsystems, described by a random unitary matrix of size N b , distributed 
according to the Haar measure. In this way, independent random unitary matrices 
describe unknown interactions in each node, which are assumed to be generic. The 
above discussion can be summarized in the following formula, describing an element 
from the ensemble: 

(2) |*> = 

where \&)fj are maximally entangled states and Uy are independent Haar unitary 
operators at each vertex. 

Our assumptions differ therefore from the model analyzed in [18], in which edges of 
the graph denote maximally entangled states of two qubits, while the vertices represent 
deterministic local unitary gates or local measurements. A more general graph model 
of quantum networks was investigated in [ij[T3] and later studied in context of entangle- 
ment percolation |19J . In this version of the model any edge represents a given bipartite 
state of two qubits, while a vertex denote a deterministic unitary swap gate or a local 
measurement. Note that in the system investigated here each dot at the end of an edge 
of a graph represents a quantum subsystem described in N dimensional Hilbert space. 
The parameter N is arbitrary, but our results are obtained under the assumption that 
the dimension is large. 
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Let us divide the set of n = 2m particles into two disjoint subsets, labeled by S and 
T. The subset T will correspond to subsystems over which the averaging is performed: 
a marginal of the graph state |^ f )(^ r | defined by the partial trace over these subsystems, 

(3) p s = Tr T \^)(% 

forms a mixed state ps supported on the remaining subspaces S = {l,...,n} \ T. 
Note that any graph T and its partition "Ptrace = {S, T} determines an ensemble of 
random mixed states ps- In [17] we studied statistical properties of such an ensem- 
ble, which depends on the topology of the graph and on its partition. Ensemble 
of random mixed states can be characterized by the average von Neumann entropy 
KH(ps) = — ETr ps log ps- By definition this quantity is equal to the average entropy 
of entanglement of the random pure state with respect to the prescribed partition 

"Ptrace = \S, T} . 

As a first example for this graphical notation, consider the graph shown in Fig. 
1(a), which consists of m — 10 edges (loops and multi-edges are allowed) and k = 5 
vertices. Fig. 1(b) shows all n = 2m = 20 subsystems denoted by small black dots and 
the random interaction between some of them represented by the gray circles at the 
vertices of the graph. The partition {S, T} is indicated graphically by placing diagonal 
crosses on the small black dots corresponding to elements in T ("traced subsystems"), 
while elements in S have no crosses on top of them ( "surviving subsystems" ) . 
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Figure 1. An example of a graph state (a) and one of its marginals (b). 



Let us now consider another example, the graph consisting of a single edge, m = 1 
presented in Fig. [2J The only edge forms a loop here, so there is a single vertex V\ only, 
k = 1. The partition splits the n = 2 partite system into two subsystems S and T. In 
this case we will relax for a moment the assumption that the sizes of both subsystems 
are equal and will denote them by iV and N', respectively. In such a case the only 
parameter of the model is the ratio between the sizes of both subsystems c = N/N'. 

As both subsystems are coupled at the vertex V\ by a random unitary matrix U G 
U(NN') the assumption on the existence of the maximally entangled state becomes 
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(a) (b) 



Figure 2. A single loop graph in the standard notation (a) and simpli- 
fied notation (b). This degenerate graph represents bipartite system (S 
- the black dot and T - crossed dot) and leads to random mixed states 
with level density described by Marchenko-Pastur distributions. 

irrelevant here. The partial trace over the subsystem T leads to a random mixed state ps 
of size N [2D] . It can be represented by a normalized Wishart matrix p$ = GG* /TrGG\ 
where G is a rectangular non-hermitian random Ginibre matrix of size N x N' (recall 
that a Ginibre matrix has i.i.d. complex Gaussian entries). The spectral density of the 
mixed state ps of size N is known to be described asymptotically by the Marchenko- 
Pastur distribution ir c (x), in the following sense 



(4) almost surely, lim pj^ = n c , 

TV— >oo 

where p^ is the empirical distribution of the rescaled eigenvalues of ps 

N 

( 5 ) UN = /^cNXjjps), 

i=l 

and the parameter c is the asymptotic ratio of dimensions c = liniAr^oo N'/N. The 
Marchenko-Pastur distribution is given by 



^4c - (x - 1 



c) 2 



(6) vr c = max(l - c, 0)<5 H — 1 [1+^2^,1+^+2^] {x) dx. 

The average von Neumann entropy of a random mixed state ps is given by the integral, 

(7) EH(ps) = kiN - J xlogx dvr c (x) + o(l) = \nN - h c + o(l), 
where h c is the entropic correction 



h 



c 



~ + clogc if c ^ 1; 

if < c < 1. 



Depending on the desired amount of generality, we are sometimes going to work on 
the model in which all subsystems are described by Hilbert spaces of the same size N. 
We will also consider the general version of the model p2] , in which only the pairs of 
subsystems connected by an edge, which describes a maximally entangled state have the 
same dimensions. Hilbert spaces of different dimensions, as in the Marchenko-Pastur 
case treated above will be allowed. However, we shall ask that these dimensions grow 
at fixed ratios, imposing the asymptotic regime dim"Hj = diN, for some fixed positive 
constants dj. 
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3. Exact area law for adapted partitions 

In this section we show that the area law holds exactly for graph states, provided that 
the marginal under consideration satisfies a particular condition, called adaptability. 

Recall that to any graph state we associate two partitions of the set of n = 2m 
subspaces: a vertex partition V ver t e x which encodes the vertices of the graph, and a 
pair partition "P e d gc which encodes the edges (corresponding to maximally entangled 
states). More precisely, two subsystems %i and %j belong to the same block of "Pvcrtcx 
if they are attached to the same vertex of the initial graph. Each edge (i, j) of the graph 
contributes a block of size two {i, j} to the edge partition P e d g e- Recall that a marginal 
(J2J) of a random graph state |^ r )(^ r | is specified by a 2-set partition P t race = {S, T}. 

Definition 3.1. A marginal p$ is called adapted if 

(9) V trace 7* 'vertex 

for the usual refinement order on partitions. In other words, a marginal is adapted if 
and only if the number of traced out systems in each vertex is either zero or maximal. 
If this is the case, then the partition boundary, which splits the graph into parts {S, T}, 
does not cross any vertices of the graph. 

Because of the above property, for adapted marginals, we can speak about traced out 
vertices, because if one subsystem of a vertex is traced out, then all the other systems 



of that vertex are also traced out. For the graph state in Figure 1(a) we consider the 
adapted marginal obtained by partial tracing vertices V 2 and V4, see Figure 3(a) 

We now define precisely what we mean by area laws [3J. The partition {S, T} defines 
a boundary between the set of vertices that are traced out and vertices that survive. In 
the subsequent sections, the following definition will be generalized to take into account 
non-adapted marginals. 

Definition 3.2. The boundary of the adapted partition {S, T} is defined as the set of 
all (unoriented) edges e = {is,jr} in the graph state with the property that is G S and 

jx G T . Equivalently, it is the number of edges of the type » •. The boundary 

of a partition shall be denoted by dS . 

The area of this boundary is its cardinality \dS\, i.e. the number of edges between S 
and T . 

In the example of Figure [31 the boundary is represented by a dashed (green) line. 



The area of the boundary in this case is 5: in Figure 3(a) , the boundary line intersects 
5 edges. 

The main result of this section is that the area law holds exactly for adapted marginals 
of graph states, where we allow arbitrary dimensions of subsystem. 
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Theorem 3.3. Let ps be an adapted marginal of a graph state \^)(^\ 
entropy of ps has the following exact, deterministic value: 



Then, the 



(10) 



Hips) - log | J] d ' N = |asi log # + log J] du 

,{i,j}£ds J \{i,j}eas 

for each value of the size parameter N. In the particular case where all the Hilbert 
spaces have dimension N (i.e. di = 1 for all i), the area law takes the form 

(11) H( Ps ) = \dS\\ogN. 

Proof. At fixed N, the random density matrix p$ has the following simple expression 



(see Figure 3(b) for a diagrammatic picture in a particular situation) 
(12) 



Ps 



Uc 



i i.it \ i, 



edge 



diN 



edge 



(g) U c 



where n^ ertex is the set of surviving vertices, Ii^ ge is the set of edges connecting surviving 
vertices and 11°^ = OS is the set of edges connecting surviving with traced vertices, 
i.e. crossed edges. Since the entropy is unitarily invariant and additive with respect to 
tensor products, it is immediate that 



(13) 



H{p s )= M<W 

{i,j}&dS 



| dS | log iV + log 




□ 



For the system corresponding to the graph shown in Figure 3(a) with all subsystems 
of size N the von Neumann entropy reads 

(14) H{p s ) = 51ogiV. 

This follows from the fact that ps is in this case a unitary conjugation of a maximally 
mixed state of size iV 5 with an arbitrary pure state of size N 6 . A graphical represen- 
tation of the reduced density operator after simplifying the unitary matrices which are 



traced out is shown in Figure 3(b) 



4. One-vertex marginals 

We shall look now at the simplest situation for which an approximate area law holds. 
We are considering marginals with a unique surviving vertex which may contain traced 
out subsystems. We refer to the results in [T7], Section 6.3, where the exact same situa- 
tion was studied. In this particular setting, let us introduce some appropriate notation 
for the subsystems of the surviving vertex V. The subsystems of V are partitioned, on 
one hand, into surviving S subsystems and traced-out subsystems T' (note that T' is 
a subset of the set of all traced systems in the original graph). Moreover, the edges 
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Figure 3. An adapted marginal for the graph state in Figure 1(a) and 
the corresponding density matrix p$. The dashed (green) line represents 
the boundary between the traced-out subsystems T and the surviving 
subsystems S. 



of the graph introduce a different partition of V, into subsystems attached by loops 
of V (denoted by F) and subsystems connected to other vertices, which form a set G. 
Hence, three important parameters with respect to the only surviving vertex: \G\, the 
number of edges connecting this vertex to other traced out vertices, \T'\ the number of 
traced out Hilbert spaces in this vertex and \S\, the number of surviving Hilbert spaces. 
For the example presented in Figure SJ one has G = {1,6,7,8}, T' = {1,2,7,8} and 
S = {3,4,5,6}. 

Except for some trivial, degenerate situations (S = V, T = V or G = V), the entropy 
of the reduced density matrix is a genuinely random, and its asymptotic behavior can 
be inferred from Theorem 6.4 of |17j . 

Theorem 4.1. In the limit N — > oo, the average von Neumann entropy of the reduced 
density matrix p$ has the following behavior: 

(\og(d s N\ s \) if\S\<\T'\ + \G\; 

(15) EH(p s ) = o(l) + <^ log(d T ^ G iVl r 'l+l G l) - h dTldG/ds if \S\ = \V\ + |G|; 

\\og(d T ,d G N\ T '\ + \ G \) if \S\ > \T'\ + |G|. 
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Figure 4. Example of a graph state marginal with only one surviving 
vertex, both in the simplified and usual notation. 



where h c is the entropic correction for a Marchenko-Pastur distribution, defined in 
equation (JHJ). In the particular case where di = 1 for all i, one has 

(16) EH(p s ) = mm{\S\, \T'\ + |G|} log N - ±6 mTWG \ + o(l). 

At this point, it is not obvious how to give an interpretation of Theorem 14.11 and of 
equation ( JTBl in particular, as an area law. Recall that because of the random unitary 
matrix acting on the vertex V, the exact indices of the traced subsystems inside the 
vertex V are irrelevant (at least in the case of trivial relative dimensions). Hence, the 
topological notion of boundary separating the sets S and T becomes ambiguous. In the 
following section, a more general definition of the boundary of a partition will be given 
via a combinatorial optimization problem. 



5. Defining the boundary surface for a general partition 

Before stating and proving the area law in the most general setting, we have to 
introduce a well-defined notion of boundary surface for a general, possibly non-adapted, 
partition {S, T}. In this section, we shall restrict our attention to the case where all 
the subsystems have the same dimension N. We shall come back to the general case of 
arbitrary ratios in Section |HJ 

A naive attempt at a definition is to consider the topological notion of boundary 
that was used in Section [3j In Figure two a priori different marginals of the same 
graph state are represented. Given the fact that independent random unitary matrices 
act on vertices V\ and V 2 , one can swap the two Hilbert spaces (say of V\) and leave 
the distribution of the reduced state invariant. Hence, the two states should have the 
same statistical properties, although they have different boundary areas: note that the 
boundary line does not cross any edges on the left hand side picture (null boundary), 
whereas the right hand side picture has a boundary area of 2 (two intersections). Hence, 
it is obvious that such a definition of the boundary area is not suitable. The rest of this 
section is devoted to introducing a new definition for the boundary area of a partition 
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motivated by a strong connection to the maximal flow problem. Note that we restrict 
ourselves to the case of identical systems, di — 1, so that dim'Hi = N. 



(a) (b) 

Figure 5. Different assignment of "crosses" yield different boundaries. 
The partition on the left has zero boundary (no cuts) and the partition 
on the right has a boundary of 2: two edges of the graph are cut by the 
dashed (green) lines. 

Consider, as before, an unoriented graph T = (V, E) having possibly multiple edges 
and loops. Consider also a counting function s : V — > N, such that, for all vertex v £ V, 
s(v) ^ deg(w) and put t(v) = deg(f) — s(v). These functions will represent the 
number of surviving (resp. traced) subsystems inside a given vertex. To the pair (r, s) 
we shall associate two combinatorial structures: a network Nr,s and a set of marked 
fattened graphs Jt, s - In Theorem 15.21 we show that the maximal flow in the network 
Nr, s is equal to maximal number of crossings in J-r>j providing the connection between 
statistical properties of the reduced state ps and a combinatorial object related to the 
partition {S, T}. This will be the main ingredient in the proof of the general formulation 
of the area law, Theorem 16.11 

Let us first describe the network Ar, s - Start with the natural network associated to 
T, with the same vertex set as T and with capacities C given by the formula 



(17) C(v,w) 



{number of edges between v and w in T if v ^ w; 
if v = w. 

To this network, add two distinguished vertices, which were called in [17] id and 7. The 
remaining capacities are defined by C(id, v) = t(v) = deg(v) — s(v) and C(v, 7) = s(v). 
The network Ar, s defined in this way has been shown in (TTJ Section 5.2] to be intimately 
connected to the statistical properties of random graph states. 

We move now to the second combinatorial object associated to the pair (r, s). First, 
starting from T, construct the fattened graph Tf at with vertex set 

(18) Vf a t = |J {vi, v 2 , ■ ■ ■ , v deg{v) } 

and with edge set Ef at corresponding to the edges of T in such a way that every two 
edges are now disjoint. We keep track of the fattening operation by a projection map 
/ : Vfat — ^ V defined by f{vi) = v. The fattening operation is depicted in Figure [6) 
A marking of fattened graph Tf at compatible with the counting function s is a subset 
M C Vf a t such that, for all v G V, 

(19) \Mnf~ 1 (v)\ = s(v). 
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We write M s to denote the fact that the marking M is compatible with the counting 
function s. Note that there always exist compatible markings and that the marking 
is unique if and only if s(v) is either or maximal for every vertex v. This extremal 
situation corresponds to adapted marginals, which were studied in Section El The 
number of crossings of a marking M, denoted cr(M), is the number of edges in E{ &t 
which connect a marked vertex with an unmarked one 

(20) cr(r fat ,M) = \{i,j} G EM | (z G M,j t M) or (z £ M,j G M)}|. 

For the first marking in Fig. Eb, one has cr(Ff at , Mi) = 6, while for the one shown in 
Fig. [6b, the number of crossings reads cr(r fat , M 2 ) = 8. 




Figure 6. The fattening of the graph in Figured] (a) and two markings, 
Mi in panel (b) and M 2 in (c), compatible with the count function (num- 
ber of marked subsystems in each vertex), s(Vi) = s(Vz) = s(Vs) = 1, 
s(Yi) = 3, and s(V&) = 2. Marked subsystems are represented by empty 
dots, and the number of crossings cr is equal to the number of edges with 
one vertex filled and the other one empty. 
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The main result of this section is the graph-theoretical Theorem 15. 2\ which makes 
use of the following important definition. 

Definition 5.1. For a graph F and a partition {S, T} of its associated nodes, define 
the area of the boundary (or, simply, areaj of the partition {S, T} as 

(21) \dS\ = maxcr(r /ai ,M), 

where the maximum is taken over all the markings M of the set of vertices of the fattened 
graph Tf at compatible with the counting function s. 

Theorem 5.2. For any graph T, the maximal flow Xr,s ^ n the network Nr,s is equal 
to the area of the boundary of the partition {S, T} 

(22) X r , s =\dS\. 

Proof. We shall prove inequalities in both directions. First, consider a compatible 
marking M s. We shall construct a set of augmenting paths in the network having 
a total flow of cr(r fat , M). For every crossing edge attached to a single vertex v of T, 
consider the augmenting path id — > v — > 7. For a crossing edge e = (v,w), where the 
black dot is in vertex v and the empty dot is in vertex w, consider the augmenting path 
id — > v — > w — > 7. In this way, to each crossing edge, we associate a unit of flow from 
id to 7, proving thus cr(Ff at ,M) ^ Xr, s . Maximizing over all compatible markings M 
proves the first inequality. 

Let us move now to proving the other direction. To this end, consider a set of 
augmenting paths in the network Ar, s , achieving the maximal flow X^ tS - Let 

id — > v 1 — )• v 2 — > ■ ■ ■ — > v k — > 7 

be such an augmenting path of length k ^ 1. If k = 1, choose an edge (vi,v{) in Ff at 
and mark one of its vertices as filled and the other one as empty. Otherwise, one can 
find the edges (t>i,f 2 ), (^2,^3), • • •, {vk-u v k) i n the fattened graph r fat . Color these 
edges in the following way: 

• (^1,^2) : vi filled, V2 empty; 

• (^2,^3) : v 2 empty, v 3 empty; 

• • ■ ■ 

• (v k -i,v k ) : v k -i empty, v k empty. 

In this way, for each augmenting path of unit flow (one can always assume this, at the 
cost of repeating edges), one assigns a unique crossing in the fattened graph. It follows 
that, for this marking M, one has cr(Ff at , M) = Xp )S . This proves the theorem. □ 

6. A GENERAL AREA LAW FOR GRAPH STATES 

This section contains the proof of the main result of the paper, Theorem 16.11 

Theorem 6.1 (Area law for random graph states). Let ps be the marginal {S, T} of a 
graphs state T. Then, as N — >■ oo ; the area law holds, in the following sense 

(23) EH( PS ) = \dS\ logiV - h r ,s + o(l), 
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where \dS\ is the area of the boundary of the partition {S, T} defined in \5.1\ and hr,s is 
a positive constant, depending on the topology of the network Ar s ( a> n d independent of 
N). 

Proof. The idea is to combine moment computations from [IT] for the random matrix 
Ps with the combinatorial identity proved in Theorem 15.21 Recall the following moment 
formula from [17J Theorem 5.5]: 

(24) Vp ^ 1, ETr(p^) = N-*™^ (\B P \ + o(l)), 

where B p C NC(p) k is the subset of non-crossing partitions (or, equivalently, geodesic 
permutations) corresponding to the augmenting paths leading to the maximal flow 
X Tt s m the network A/r,s (for details, see [H]). We can restate the above asymptotic 
expression as a limit: (we write simply X = Xr,s) 

(25) \/p > 1, jim E^Tr [(N x p s ) p ] = \B P \. 
In other words, the measures 

1 N x 

(26) p N = E— x 6 N*\ t 



N x 

i=l 



(PS) 



have limiting moments given by \B p \. Note that there is a unique probability measure 
having moments \B p \, since one has the bound \B P \ ^ Cat J ^ (4 fc ) p which is exponen- 
tial, and thus, by Carleman's condition, these moments uniquely define the probability 
measure. Since the limit points of the tight sequence of measures p^ are uniquely 
determined by their moments, it follows that [221 Theorem C.9] pn converges weakly 
towards a measure p satisfying 

Finally, one has 

(27) EH(p s ) =X\ogN- h r>s + o(l), 

with a correction term equal to the entropy of the asymptotic measure p, 

(28) hr,s — J x log xdp(x). 

□ 



If the measure p^ converges to the Marchenko-Pastur distribution with parameter 
c, the above term coincides with the entropy h c defined in (JHJ). 
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7. Rank of random graph states and a transport problem 

We start by looking at a linear algebra problem, the maximum rank of a marginal of 
a graph state, over the set of unitary operations on vertices. 

Theorem 7.1. The maximum rank of a graph state marginal ps is the area of the 
boundary of the partition {S, T} 

(29) max rkp 5 = \dS\. 

v ' u u ...,u k r 1 1 

Moreover, this maximum can be achieved by choosing the Ui to be permutation matrices 
and the state p* s which achieves the maximum can be taken maximally mixed. 

The above result shows that the maximum achievable rank of a graph state is equal, 
asymptotically, to the rank of a random marginal. We shall investigate further this 
property by presenting an alternative formulation to our problem, bridging together 
the max-flow, max-crossings and the rank aspects. It also has the advantage of being 
"operational". The notation below mirrors the one in the rest of the paper and in |17j . 
Note that there is no randomness in the problem below, and N can be arbitrary. We 
consider below the simplest one-qubit case, N = 2. 

Consider the following problem (see Figure [7] for an example). 



Fx (2,1) V a (l,3) V 3 (2,l) 

i • • • | • • • 




(c) 



Figure 7. An instance of the problem, where Si and T, are noted in 
parenthesis after each Vi. In (b), we create locally singlet states. The 
states sent to A are circled in red in (c), the other ones being sent to B. 
Note that at each site we consider Si + T particles. 

Problem. The company EntanglementFactory has k research facilities around the 
world, let us call them Vx, . . . , T4- As a result of past experiments, entangled states 
are shared between pairs of these laboratories, as follows : labs Vi and Vj share 
singlet states |$ + ) G C^CgxC^. Each facility Vi has an unlimited supply of extra iV-dits 
which are not entangled with anything else. One day, company EntanglementFactory 
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receives an order from the company WantEntanglement for a supply of entangled states. 
Company WantEntanglement has two research facilities A and B and offers to ship Au- 
dits from each factory V* to A and/or B, as follows: Si Audits can be shipped from Vi 
to A and Tj Audits can be shipped from Vi to B. Company WantEntanglement pays 
$10 for each unit of entanglement it will have between stations A and B. What is 
the maximal profit the company EntanglementFactory can make, just by using local 
unitary matrices at each site Vi ? 

We shall answer this question in the following three scenarios, each situation imposing 
some physical restrictions or liberties on the system. 

Scenario 1 : No initial entanglement. Suppose that all the entangled particles 
shared between pairs of VJ's are lost. The best we can do is to create locally, at each 
Vi, maximally entangled (singlet) states and to ship one half to A and the other half to 
B. The number of entangled pairs between A and B will then be 

Yi = ^2mm(Si,Ti). 

i 

Scenario 2 : Global operations are allowed. Suppose that company Entangle- 
mentFactory has the ability of performing global operations on all of its facilities Vi. 
Then it is easy to produce a maximally entangled state between A and B: 

Y 2 = min(J] S^Ti). 

i i 

Scenario 3 : Entanglement with local operations. Without any assumptions, we 
show that 

First, notice that the problem in this case is just as a restatement of the rank theorem 
discussed earlier, before taking the partial trace of the state. The facilities A and B 
define the partition with respect to which the partial trace is considered. The entropy 
of entanglement of the pure state shared between A and B is just the von Neumann 
entropy of the reduced density matrix p$. Note that the following inequality concerning 
the number of entangled pairs in different scenarios holds Y ± ^ Y 3 ^ Y 2 . 

One can characterize entanglement by using the generalized Renyi entropy H q (p) : = 
lnTrp 9 . Note that in the limit the Renyi parameter q tends to unity this expression 
reduces to the von Neumann entropy, lim,^! H q (p) = H(p). Furthermore, the rank rk 
of a matrix is given by the generalized entropy of order zero, logrkp = H (p). 

To summarize, the results on the rank and on the Renyi entropy H q have the following 
translation 

Theorem 7.2. There exist local unitary operations Ui (which can be taken to be per- 
mutation matrices, i.e. at each site Vi we just have to say where each particle goes, to 
A or B) such that 

H (p) = H^p) = H p (p) = X r , s logiV, 
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for any q > 0. In other words, p is essentially a maximally mixed state. Moreover, the 
permutation matrices involved can be computed efficiently, using a flow algorithm. 

In the random case (suppose the engineers of company EntanglementFactory are 
on vacation, so the staff decides to implement at each site random, independent 
local unitary transformations), we have the following result, which is a restatement of 
Theorem 16.11 

Theorem 7.3. For independent random Haar unitary matrices Ui, we have, almost 
surely as N — > oo ; 

H (p) = X r , s \ogN 

and 

EH^p) = X r , 5 logiV - h r ,s + o(l). 

Comparing the two results above, one concludes that the random choice is nearly 
optimal. 



8. Some results for different subsystem dimensions 

In this final section, we analyze simple graphs in the general setting, where we allow 
subsystems to have different dimensions. Although we can not state a general area 
law for marginals of such graph states, we perform direct computations in some simple 
settings using the full machinery developed in [I7j . Our main tool is the following 
result, valid for random graph states with subsystems of dimension diN . 

Theorem 8.1 (see [T71 Theorem 5.4]). The asymptotic moments of a graph state mar- 
ginal ps are given by the formula 

k k 

(so) ETr( P p s ) = (i + (i))j\r*<^ £ n (^) #(7 " ft) n 

(^l,.,fe)£B»=l »=1 
l<i<j<fc »=1 

where B is a set of permutations that can be computed from the network Afr,s associated 
to the marginal and X is the maximal flow id — > 7 in the same network. For each block 
Ci of Yl vertex , we write Si = S D Ci and Tj = T fl Cj. Finally, Eij is the set of edges 
going from vertex i to vertex j . 

The above general expression for the average p-th moment allows us to analyze ensem- 
bles of quantum states corresponding to the following exemplary graphs, each containing 
4 subsystems, but having different geometry. 
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Black holes graph. We are going to discuss here a simple graph described by two 
edges joined in one vertex, used to model the trans-horizon entanglement during the 
process of evaporation of a black hole (23j[23] . The entire system is thus composed 
out of four subsystems, two of which have the same dimension equal to d\N, while the 
other two have the dimensions (I2N. The ratios are treated as parameters of the 
model. We shall consider two marginals of this graph state, see Figure [HJ In both cases, 
the network associated to the marginal is the same, and has a maximum flow X = 2. 
Moreover, the set B of permutations achieving this maximum flow is 



(31) 



B = {(ft, A,, AO e S 3 P : id = ft ^ ft < & = 7}- 



diiV d 2 N 

o- -o- -O 
Vi v 2 v 3 

(a) 



v 2 

(c) 





Vi 




r 2 

(b) 




Va 




Vi 



r 2 

(d) 



Va 




(e) 



Figure 8. The black hole graph, in the simplified notation (a) and the 
diagram with 4 subsystems (b). In the middle row, two of its marginals: 
on the left (c), G?iiV-sized subsystems are traced out and on the right (d), 
subsystem of both dimensions are traced out. In the bottom row, the 
network associated to both marginals (e). 
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First case : subsystems of size diN are traced. Applying Theorem 18.11 to this setting, 
we obtain 

(32) ETr p p s = (1 + o(l))N-^ £ df^df^df (ft) rff m 

id=/3isC^ 2 s:/33=7 

.df^-'dt^^-^dldiyv 



a£NC(p) ^ 2 ' 



1 / rl \ 2#a 

lim -J— ETr( dlN 2 Ps ) p = y (- 

N^ood 2 N 2 V 1 Hb) ^ \d 2 ) 



In other words 
(33) 

\ — ■ X (17-,, \ - 

<reJVC(p) 

Thus, the rescaled random matrix d\N 2 ps G ./^^(C) converges in moments to the 
Marchenko-Pastur distribution of parameter d\jd\. Let us now compute the average 
entropy of this random matrix. Use 

(34) jim ^^H{d\N 2 ps) = -h^ 
to show that 

f iog(^iv 2 ) -A iid 1 ^d 2 

(35) EH(p s ) = (1) + 9 ? 

\log(d 2 iV 2 )-|| iid 1 <d 2 , 

= log(rf 2 iV 2 )-^ + o(l), 

where d = min((i 1 , d 2 ) and -D = max(c?x> ^2)- Note that the formula above is symmetric 
in d\ and d 2) a consequence of the fact that the non-zero spectra of the two reduced 
density operators of a pure state are identical. 

Second case : subsystems of both sizes are traced. In this case, the moment formula 
reads 

(36) ETrp£ = (l + o(l)) N- 2 ^- l \d 1 d 2 y p+1 J2 L 

id=/3 l5 :/32^/33=7 

Thus, the rescaled random matrix d\d 2 N 2 p G M-d^N 2 ^) converges in moments to the 
Marchenko-Pastur distribution of parameter 1. The entropy computation in this case 
is easier: 

(37) EH(ps) = log(cM 2 iV 2 ) - i + o(l). 
Note that the two entropy formulas agree in the case d\ = d 2 . 
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Double line graph - the oxygen molecule O2. The graph presented in Figured], 
which can be symbolically represented by O = O, might be interpreted as the oxygen 
molecule. We will discuss here the general version of the model in which there are 
two pairs of subsystems of size d\N and d%N respectively, and look at two different 
marginals. The two marginals correspond to the same network, which has a maximal 
flow X = 2. The set B of permutations achieving this maximum flow is 

(38) B = {(/3i,/3 2 ) G S 2 p : id < (3 1 < < 7 and id < fa < 0i < 7} 

= {(/?!, (3 2 )£S 2 p : id^/3 1= /3 2 ^ 7 }, 

which is in bijection with the set found for the "black-hole" graph. 




FIGURE 9. The oxygen graph, in the simplified notation (a) and the 
diagram with 4 subsystems (b). In the middle row, two of its marginals: 
on the left (c), c?iiV-sized subsystems are traced out and on the right (d), 
subsystem of both dimensions are traced out. In the bottom row, the 
network associated to both marginals (e). 



First case : subsystems of size diN are traced. Proceeding as in the case of the "black- 
hole" graph, we obtain the same result as in the first case above. Thus, as before 

(39) EH(p s ) = log(d 2 iV 2 ) - 5^5 + 
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where d = min(d 1 , d 2 ) and D = max(d 1 , d 2 ). 

Second case : subsystems of both sizes are traced. Again, we obtain the same result as 
in the corresponding case of the "black-hole" graph 

(40) EH(p s ) = \og(d 1 d 2 N 2 ) - X - + o(l). 

To conclude, note that in general, the relative size of the traced out systems matters 
(unless di = d 2 ). At this point, we have no interpretation for the surprising fact that 
the two graphs studied here yield the same output entropies. 
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